Generic features of the fluctuation dissipation relation in coarsening systems 
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The integrated response function in phase-ordering systems with scalar, vector, conserved and 
non conserved order parameter is studied at various space dimensionalities. Assuming scaling of 
the aging contribution Xaaiijiw) ~ we obtain, by numerical simulations and analytical 

arguments, the phenomenological formula describing the dimensionality dependence of in all 
cases considered. The primary result is that vanishes continuously as d approaches the lower 
critical dimensionality di . This implies that i) the existence of a non trivial fluctuation dissipation 
relation and ii) the failure of the connection between statics and dynamics are generic features of 
phase ordering at dz,. 
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After the groundbreaking work of Cugliandolo and 
Kurchan P] on mean field spin glasses, the study of the 
out of equilibrium linear response function has been gain- 
ing an increasingly important role in the understand- 
ing of slow relaxation phenomena. The key concept is 
that of the fluctuation dissipation relation (FDR) 2] . In 
terms of the integrated response function x(<,tu,), i.e. 
the response to an external field acting in the time in- 
terval {tw,t), an FDR arises if x(*j^i«) depends on time 
only through the unperturbed autocorrelation function 
C{t,tw). If this happens, there remains defined a func- 
tion X — S{C) which generalizes the fluctuation dissipa- 
tion theorem into the out of equilibrium regime. 

The existence of an FDR is important for several rea- 
sons 0. Here we focus on a specific aspect: to what 
extent the FDR shape is revealing of the mechanism 
of relaxation and of the structure of the equilibrium 
state. In particular, we aim at dispelling the common 
belief that relaxation by coarsening and a simple equi- 
librium state do necessarily imply a fiat or trivial FDR, 
i.e. S{C) = 1 — qsA when C falls below the Edwards- 
Anderson order parameter qEA- 

To appreciate the relevance of the problem, consider 
that, by reversing the argument, the observation of a 
non flat FDR, would rule out relaxation by coarsening. 
This is a statement of far reaching consequences. For 
instance, an argument of this type plays a role in the 
discrimination between the mean field and the droplet 
picture of the low temperature phase of finite dimensional 
spin glasses • In that case, the flnal conclusion may well 
be right, but for the argument to be sound, the behavior 
of the response function, when relaxation proceeds by 
coarsening, needs to be thoroughly understood. 

As a contribution in this direction, we have undertaken 
a large program of systematic investigation of the FDR in 



the context of phase ordering systems jj] , which provide 
the workbench for the study of all aspects of relaxation 
driven by coarsening. We have considered pure ferromag- 
netic systems quenched from above to below the critical 
point. We have covered the whole spectrum of systems 
with non conserved (NCOP), conserved (COP), scalar 
[N = 1) and vector [N > 1) order parameter at differ- 
ent space dimensionalities d, where N is the number of 
components of the order parameter. The manifold of the 
systems considered is displayed in Table I. Some of these 
(marked by a dot) have been studied before. The im- 
portant novelty is that, with the new entries, the picture 
becomes rich enough to promote to generic the behav- 
ior previously observed in the particular case of the Ising 
model 0, S 0i 13 and in the large N model • Namely, 
that a flat FDR is obtained for d > d-L, while a non flat 
FDR is found for d = dL, where is the lower critical 
dimensionality. The implication is that a flat FDR is not 
a necessary condition for coarsening. 

To explain in more detail, let us recall that, quite 
generally, one can write x{t, tw) = Xst{t-t^)+Xag{t, t^)- 
The first is the stationary contribution due to the fast de- 
grees of freedom which rapidly equilibrate with the bath, 
while the second is the aging contribution coming from 
the slow out of equilibrium degrees of freedom. What one 
can also show, in general, is that a fiat FDR is obtained 
if Xagit, tyj) vanishes asymptotically Now, in phase 

ordering for large tw one expects the scaling behavior 



(1) 



from which it follows that the FDR is or is not fiat ac- 
cording to > or flp^, < (here we will restrict to 
the case with = 0). Therefore, investigating the FDR 
shape requires the investigation of a^. 

Let us see what is the situation with this exponent. In 
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the Ising model 0, 0, 0, II] and in the large N model , 
we have found that depends on dimensionality accord- 
ing to 
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where S enters the time dependence of the density of 
defects and di, du are the two special dimensionalities 
(with du < du) where — and above which = 
6, respectively. The density of defects goes like p{t) ^ 
L{ty ~ where L{t) - i^/^ is the typical defect 
distance, z is the dynamic exponent and n = 1 or n = 2 
for scalar or vector order parameter ^j]. Hence, 6 = n/z. 

Opposite to Eq. |2l stands a qualitative argument ac- 
cording to which Xag{t,tw) ought to be simply propor- 
tional to the defect density 0, 0|. Namely, = 5 
at any dimensionality. There are no measurements or 
derivations of supporting such a statement. Further- 
more, the arg ument is incompatible with the exact result 
Ox ~ 0, ^3 i'^ the d — 1 Ising model. The implication 
for the FDR is that it should be always flat, except for 
possible pathological cases. This point of view is diame- 
trally opposed to the one underlying Eq. (0), where the 
vanishing of at d = di, far from being a pathology, 
comes out as a limiting behavior in the smooth dimen- 
sionality dependence of for d < du- 

The understanding of the FDR in coarsening systems 
requires to clarify whether a^^ does or does not to depend 
on d. Here we present strong evidence supporting Eq. (O 
as the generic pattern of behavior. 

We have computed x(i, tw) for systems quenched from 
infinite to zero final temperature. This is computation- 
ally efficient and can be done without loss of general- 
ity, since all quenches below Tq are controlled by the 
T = fixed point 0. In all cases we have used the 
time dependent Ginzburg-Landau equation Q, except 
for NCOP with iV > 1 and d > 2 where the Bray- 
Humayun algorithm has been used 0|. The sta- 
tionary response has been computed from equilibrium 
simulations. The aging part, then, has been obtained 
from Xag{t,tw) = x{tiiw) ~ Xst{t - t^)- To get a^, 
one ought to extract the dependence of Xag{t,tw) for 
fixed X = t/tuj However, this is computationally 

very demanding and would make it impossible to get the 
vast overview we are aiming at. So, we have measured 
from the large t behavior for a fixed t^, assuming 
Xag{t,tw) ~ t~"->^. This holds if x{x) ~ x~°-^ for x > 1, 
which has been verified in the NCOP scalar case 0, Q 
and it is an exact result in the soluble models 0, Q . The 
assumption is that it holds in general. The choice of 
is inessential provided it is larger than some microscopic 
time necessary for scaling to set in 0. 

The time dependence of Xag{t,tw) is depicted in 
Figs, n 121 and 13 We have extracted from the asymp- 



totic power law decay and we have collected all results, 
old and new, in Table II. At du we have used the para- 
metric plot Xag{C) (insets of Figs, d El and EJ, showing 
more effectively the absence of asymptotic decay, due to 
— 0. In Table II, we have also reported the values 
of predicted by Eq. The comparison with the 

computed values is quite good. For convenience, we have 
collected in Table HI the values of all the parameters 
entering Eq. (0). Finally, Fig. ^ provides the pictorial 
representation of Table II, and it is the main result in 
the paper. 

Let us now comment the results. From Fig. 0] it is ev- 
ident that the pattern of behavior predicted by Eq. 
is obeyed with good accuracy in the scalar cases, with 
du — 1 and du — 3. In the vector cases, given the great 
numerical effort needed, values of N were chosen accord- 
ing to the criterion of the best numerical efficiency, to- 
gether with the requirement to simulate both systems 
with {N < d) and without {N > d) stable topological 
defects. The overall behavior of the data in Fig.^lshows 
that Eq. ^ well represents the dimensionality depen- 
dence of also in the vector case with du = 2 and 
du = 4. Finally, the insets in Figs [3 and |3 (together 
with the analogous figures for the d — 1 Ising model in 
Ref.s 0, and in the large N model 0) show quite 
clearly that = and a non flat FDR are common 
features in phase ordering kinetics at d^. 

At this stage Eq. is a phenomenological formula. 
Apart from the exact solution of the large N model Q, 
there is no derivation of Eq. Here we propose an 

argument for the dependence of on d in the scalar 
case. It is based on two simple physical ingredients: 
a) the aging response is given by the density of defects 
p{t) times the response of a single defect Xag{t,tw) = 
p{t)Xag{'t'^tw) and b) each defect responds to the per- 
turbation by optimizing its position with respect to the 
external field in a quasi-equilibrium way. In d = 1 this 
occurs via a displacement of the defect. In higher dimen- 
sions, since defects are spatially extended, the response 
is produced by a deformation of the defect shape. 

We develop the argument for a 2-d system, the ex- 
tension to arbitrary d being straightforward. A defect 
is a sharp interface separating two domains of opposite 
magnetization. In order to analyse Xagi^^^w) we con- 
sider configurations with a single defect as depicted in 
Fig. The corresponding integrated response function 
reads 6] xlgit.t^) = I / {h^ C"-^) ! dxdy {Six,y))h{x,y), 
where S{x, y) is the order parameter field which sat- 
urates to ±1 in the bulk of domains. h{x,y) is the 
external random field with expectations h{x,y) — 0, 
h{x,y)h{x' ,y') = h?5[x — x')5[y — y') and C is the 
linear system size. The overbar and angular brack- 
ets denote averages over the random field and ther- 
mal histories, respectively. With an interface of shape 
Zw{y) at time (Fig. Ell, we can write X^git^t^) = 
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-l/{h^Cd-^)J^^^ Eh Ph{{z{y)},t), where Ph{{z(y)},t) is 
the probability that an interface profile {z{y)} occurs 
at time t and Eh = - Jq dy J^^^^^^ dxh{x,y)sign[z{y) - 
^wiv)] is the magnetic energy. We now intro- 
duce assumption b) making the ansatz for the cor- 
rection to the unperturbed probability Po{{z{y)},t) 
in the form of a Boltzmann factor Pfi{{z{y)},t) = 
Po{{ziy)},t)eM-l3Eh) - Po({z(y)}, ^[1 - Then 
Xig{t,t^) = -l/ih^C''-')J^,yEh{l-f3Eh)Poi{ziy)},t). 
Taking into account that the term linear in Eh van- 
ishes by symmetry and neglecting Zyj{y) with respect 
to z{y) for t 3> tu), we eventually find fi~^Xag{t^tw) = 
'^'"''/{.}/o^^2/l^(y)l^o({^(y)},t)- This defines a length 
which scales as the roughness of the interface given 
hy W{t) = [/:i-'^/^^Jdyz(j/)2Po({z(y)},t)]i/2. The 

behavior of W{t) in the coarsening process can be in- 
ferred from an argument due to Villain [T^ . In the case 
d < 3, when interfaces are rough for NCOP one 
has W{t) - t(3-rf)/4, while for COP W{t) - t^^-d)/&^ 
with logarithmic corrections in both cases for d = 3. For 
(i > 3 interfaces are flat and W{t) ~ const. Finally, mul- 
tiplying xkg by p{t) ^ L{t)'^ Eq. (0) is recovered 1^ 
and du is identified with the roughening dimensionality 
dn = 3. The relevance of roughening in the large time be- 
havior of t^) has been independentl y p ointed out by 
Henkel, Paessens and Pleimling in Ref. |l9j- The crucial 
difference with these authors is that they believe rough- 
ening to be unrelated to aging behavior, while we claim 
the opposite. 

In summary, we have investigated the scaling proper- 
ties of the response function over a large variety of sys- 
tems designed to bring forward the generic features when 
relaxation is driven by coarsening. The primary result is 
that the exponent depends on dimensionality and that 
it vanishes smoothly as d c^l . This implies that a non 
trivial FDR is not exceptional, rather is the rule for coars- 
ening systems at dL- Another important consequence is 
that the failure of the connection between statics and dy- 
namics at di is also a generic feature of coarsening. 
The connection between the FDR and the overlap prob- 
ability function is derived \vi\ under the assumptions of 
stochastic stability and that goes to the equilib- 

rium value as i — > oo. The latter assumption does not 
hold at dL due to the existence of a non flat FDR (insets 
of Figs 11121.^ . which makes the Hmiting value of xi^iiw) 
to rise above the equilibrium value. Obviously, the im- 
portant and, as of yet, unanswered question is why all 
this happens at d^. The scaling behavior of the response 
function reported in this Letter adds to the many already 
existing challenges posed by a theory of phase ordering 
kinetics. 
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TABLE I: The manifold of systems considered. Entries with 
dots correspond to systems studied in Refs. IMIEQIS' 
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FIG. 1: Xti9(i,im) against t — t^ for A'^ = 1 with COP. Lattice 
sizes, realizations and t^: 512^^, 41 and 30 for d — 2; 128^, 39 
and 40 for d = 3; 60^, 6 and 31 for d — A. The dashed lines 
are the slopes from Eq. ||5J. In the inset parametric plot for 
d=l from Ref. Q 
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FIG. 2: Xa.g{t,tw) against t ~ with COP. Lattice sizes, 
realizations and tu>- 96'^, 89 and 35 for d = 3 and A'^ = 5; 
SO"*, 82 and 35 for d = 4 and iV = 2. The dashed lines 
are the slopes from Eq. ||5J. In the inset parametric plot for 
d — 2,N = A. Lattice size, realizations and tw- 512^, 232, 
500). 




FIG. 3: Xag{i,'tiD) against t — with NCOP. Lattice sizes, 
realizations and tw- 180^, 1445 and 2 for d = 3 and N — 2; 
140^ 1486 and 0.3 for d = 3 and iV = 5; 40^*, 486 and 0.3 for 
d = 4 and N = 6. In the inset parametric plot for d — 2 and 
A'' = 10. Lattice size, realizations and tw'- 1024^, 22 and 20). 
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FIG. 4: Survey of results in table UTI 
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FIG. 5: 
(dashed 



Configurations with a single interface at time t, 
line) and at time t (continuous line). 



